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Description 
Technical Field 

s [0001 ] This invention relates to the art of wireless communications, and more particularly, to wireless communication 
systems using multiple antennas at the transmitter and multiple antennas at the receiver, so called multiple-input mul- 
tiple-output (MIMO) systems. 

Background of the Invention 

10 

[0002] It has previously been assumed for MIMO systems that the time dispersion between one or more of the 
transmit antennas and one or more of the receive antennas was negligible, i.e., the various paths were nonresolvable, 
as often occurs in low bandwidth systems, because the pulse width is longerthan the channel time dispersion. However, 
it has been recognized that under typical urban (TU) conditions, i.e., the conditions of the so-called "TU" model, that 
15 the time dispersion between one or more of the transmit antennas and one or more of the receive antennas is non- 
negligible. Such a non negligible time dispersion causes the various paths to interfere with each other, resulting in a 
bit error rate floor, i.e. minimum, and so the resulting bit error rate is unacceptable. 

[0003] I. Ghauri and D. Slock have shown, in "Linear Receivers forthe DS-CDMA Downlink Exploiting Orthogonality 
of Spreading Codes", 32 nd Asilomar Conference, Nov. 1 -4, 1998 pp. 650-4, that a minimum mean square error (MMSE) 
20 equalizer operating on received code division multiple access (CDMA) chips can be employed to compensate for time 
dispersion in a single transmit, single receive antenna system, thus reducing the bit error rate floor and improving 
performance. 

Summary of the Invention 

25 

[0004] We have recognized, in accordance with the principles of the invention, that in a MIMO system the bit error 
rate floor caused by time dispersion can be reduced by employing a joint equalizer for all of the respective transmit 
antenna - receive antenna pairings that are possible in the MIMO system. Advantageously, the resulting joint equali- 
zation compensates not only forthe impact of the channel on the transmit antenna - receive antenna pairings but also 
30 for the interference of the other transmit antennas on any given receive antenna. In a particular embodiment of the 
invention, the joint equalizer is a joint minimum mean square error (MMSE) equalizer, and in such an embodiment the 
joint equalization outperforms simply replicating the prior art minimum mean square error (MMSE) equalizer for each 
transmit antenna - receive antenna pairing. 

[0005] In one embodiment of the invention, which is especially useful for CDMA, after the equalization is completed 
35 the resulting chip streams, one for each transmit antenna, are despread in the conventional manner and then the 
resulting dispread symbols may be further processed in the conventional manner. Alternatively, instead of further 
processing the despread symbols in the conventional manner, the despread symbols may be processed, in accordance 
with an aspect of the invention, so as to have their soft bits computed through the use of a posteriori probability (APP) 
metric. Prior to computing the soft bits, the despread symbols may be spatially whitened using a spatial whitening filter. 
40 [0006] In another embodiment of the invention, which is also especially useful for CDMA, the equalizer is iteratively 
computed so that a symbol from one transmit antenna is determined during each iteration. Initially the received samples 
are stored in a memory. After a symbol for an antenna is determined, the received samples for each of the remaining 
antennas are then recomputed by subtracting out the determined symbol from the samples as they existed prior to 
determining the symbol. Once all the symbols for all of the transmit antennas are determined for a symbol period are 
45 determined the operation begins anew with the samples corresponding to the next symbol period. 

Brief Description of the Drawing 

[0007] In the drawing: 

50 

FIG. 1 shows an embodiment a multiple-input multiple-output (MIMO) wireless system in which the bit error rate 
floor caused by time dispersion is reduced by employing a joint minimum mean square error (MMSE) equalizer 
for ail of the respective transmit antenna - receive antenna pairings that are possible in the MIMO system, in 
accordance with the principles of the invention; 
55 FIG. 2 shows an exemplary process, in flow chart form, for the overall operation of system of FIG. 1 ; 

FIG. 3 shows in more detail the process by which the weights employed by the joint equalizer of FIG. 1 are deter- 
mined; 

FIG. 4 shows another exemplary process, in flow chart form, for the overall operation of system of FIG. 1 ; 
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FIG. 5 shows another embodiment a multiple-input multiple-output (MIMO) wireless system in which the bit error 
rate floor caused by time dispersion is reduced by employing a joint minimum mean square error (MMSE) equalizer 
for all of the respective transmit antenna - receive antenna pairings that are possible in the MIMO system, in 
accordance with the principles of the invention; 
5 FIG. 6 shows a more detailed version of a buffer-subtractor of FIG. 5; 

FIG. 7 shows an exemplary process, in flow chart form, for the overall operation of system of FIG. 5; 

FIG. 8 shows a particular embodiment of the joint equalizer of FIG. 1 , in which the equalization is performed in the 

discrete frequency domain, in accordance with an aspect of the invention; and 

FIG. 9 shows a particular embodiment of the joint equalizer of FIG. 1 , in which the equalization is calculated in the 
10 discrete frequency domain and applied in the time domain. 

Detailed Description 

[0008] The following merely illustrates the principles of the invention. It will thus be appreciated that those skilled in 
15 the art will be able to devise various arrangements which, although not explicitly described or shown herein, embody 
the principles of the invention and are included within its spirit and scope. Furthermore, all examples and conditional 
language recited herein are principally intended expressly to be only for pedagogical purposes to aid the reader in 
understanding the principles of the invention and the concepts contributed by the inventor(s) to furthering the art, and 
are to be construed as being without limitation to such specifically recited examples and conditions. Moreover, all 
20 statements herein reciting principles, aspects, and embodiments of the invention, as well as specific examples thereof, 
are intended to encompass both structural and functional equivalents thereof. Additionally, it is intended that such 
equivalents include both currently known equivalents as well as equivalents developed in the future, i.e., any elements 
developed that perform the same function, regardless of structure. 

[0009] Thus, for example, it will be appreciated by those skilled in the art that any block diagrams herein represent 
25 conceptual views of illustrative circuitry embodying the principles of the invention. Similarly, it will be appreciated that 
any flow charts, flow diagrams, state transition diagrams, pseudocode, and the like represent various processes which 
may be substantially represented in computer readable medium and so executed by a computer or processor, whether 
or not such computer or processor is explicitly shown. 

[0010] The functions of the various elements shown in the FIGs., Including any functional blocks labeled as "proc- 

30 essors", may be provided through the use of dedicated hardware as well as hardware capable of executing software 
in association with appropriate software. When provided by a processor, the functions may be provided by a single 
dedicated processor, by a single shared processor, or by a plurality of individual processors, some of which may be 
shared. Moreover, explicit use of the term "processor" or "controller" should not be construed to refer exclusively to 
hardware capable of executing software, and may implicitly include, without limitation, digital signal processor (DSP) 

35 hardware, network processor, application specific integrated circuit (ASIC), field programmable gate array (FPGA), 
read-only memory (ROM) for storing software, random access memory (RAM), and non-volatile storage. Other hard- 
ware, conventional and/or custom, may also be included. Similarly, any switches shown in the FIGS, are conceptual 
only Their function may be carried out through the operation of program logic, through dedicated logic, through the 
interaction of program control and dedicated logic, or even manually, the particular technique being selectable by the 

40 implementor as more specifically understood from the context. 

[0011] In the claims hereof any element expressed as a means for performing a specified function is intended to 
encompass any way of performing that function including, for example, a) a combination of circuit elements which 
performs that function or b) software in any form, including, therefore, firmware, microcode or the like, combined with 
appropriate circuitry for executing that software to perform the function. The invention as defined by such claims resides 

45 jn the fact that the functionalities provided by the various recited means are combined and brought together in the 
manner which the claims call for. Applicant thus regards any means which can provide those functionalities as equiv- 
alent as those shown herein. 

[0012] Software modules, or simply modules which are implied to be software, may be represented herein as any 
combination of flowchart elements or other elements indicating performance of process steps and/or textual description 
50 Such modules may be executed by hardware which is expressly or implicitly shown. 

[0013] Unless otherwise explicitly specified herein, the drawings are not drawn to scale. 

[0014] In the description, identically numbered components within different ones of the FIGs. refer to the same com- 
ponents. 

[0015] FIG. 1 shows an exemplary embodiment of a multiple-input multiple-output (MIMO) wireless system in which 
55 the bit error rate floor caused by time dispersion is reduced by employing a joint minimum mean square error (MMSE) 
equalizer for all of the respective transmit antenna - receive antenna pairings that are possible in the MIMO system, 
in accordance with the principles of the invention. Shown in FIG. 1 are a) transmitter 101; b) transmit antennas 103; 
including transmit antennas 1 03-1 through 1 03-M; c) receive antennas 1 05; including receive antennas 1 05-1 through 
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105-N; d) receiver front-end 107; e) joint equalizer 109; f) optional despreaders 111 ; g) soft bit mapper 112, which may 
include optional spatial whitening filter 113; and optional a posteriori probability (APP) metric processor 115. 
[0016] Transmitter 101 is a MIMO transmitter, e.g., one in which an original data stream is divided into substreams 
and each resulting substream is transmitted as a modulated radio signal via an individual one of transmit antennas 
5 1 03. The transmitted signals pass to the receiver over a time dispersive channel such that signals from each transmit 
antenna 103 reach each of receive antennas 105. 

[001 7] Receive antennas 1 05 convert the radio signals impinging upon them into electrical signals, which are supplied 
to receiver front-end 107. Receiver front-end 107 operates conventionally to produce a stream of binary numbers 
representing samples of the radio signals received at antennas 105. Typically receiver front-end 107 performs radio 
10 frequency downconversion, filtering, sampling, and analog-to-digital conversion. The resulting samples are provided 
to joint equalizer 109. 

[0018] Joint equalizer 109 compensates for the effects of the transmit signals from each of antennas 103 having 
passed through the channel as well as the interference that results from transmitting via multiple antennas simultane- 
ously. The output of joint equalizer 109 is M, i.e., the number of transmit antennas, streams of corrected symbols, or 
15 in the case of CDMA, streams of corrected chips which when properly combined form symbols. Operation of joint 
equalizer 109 will be described more futly hereinbelow. If CDMA is employed, the output of joint equalizer 109 is. 
supplied to optional despreaders 1 1 1 , which produces symbols from the stream of chips supplied by joint equalizer 1 09. 
[0019] The symbols produced may then be further processed in the conventional manner for a MIMO system, e.g., 
soft bits may be developed by soft bit mapper 1 12 for use in a decoder, e.g., the well known "Turbo decoder 1 '. Alterna- 
te tively, the symbols may be supplied within soft bit mapper 112 to optional spatial whitening filter 113, which makes the 
noise equal on each branch. Note that the whitening is performed only in the space domain. If whitening is performed 
in the time domain some temporal dispersion will be introduced into the signal. The symbols, or whitened symbols if 
optional spatial whitening filter 1 1 3 is employed, may further be supplied to optional a posteriori probability (APP) metric 
processor 115 within soft bit mapper 112, in accordance with an aspect of the invention. APP metric processor 115 
25 performs a particular type of mapping from symbols to soft bits. Operation of APP metric processor 1 15 is described 
more fully hereinbelow. 

[0020] FIG. 2 shows an exemplary process, in flow chart form, for the overall operation of system of FIG. 1 . Prior to 
performing the process of FIG. 2, initial values of parameters M, N, L, P r E, and dmust be determined. What each of 
these parameters represents is listed in Table 1 . Additionally, prior to performing the process of FIG. 2 it is necessary 
30 to determine noise covariance R pp , which contains samples of the autocorrelation of the chip pulse shape r(t), when 
CDMA is employed, or the symbol pulse shape autocorrelation, when CDMA is not employed. 



Table 1 - Parameter Definitions 



35 


M 


number of transmit antennas 


N 


number of receive antennas 




L 


length of the channel impulse response (chips) 




P 


over-sampling factor 




E 


length of the equalizer (chips) 


40 


d 


equalizer delay (chips) 




On 2 


power of interference + noise (RF bandwidth) 




Ox 2 


power of downlink signal (RF bandwidth) 




G 


Number of CDMA chips per symbol 



[0021 ] The process of FIG. 2 is executed periodically, with a periodicity that preferably does not exceed the coherence 
time, T co , which is the time duration for which the channel properties are substantially constant. The process is entered 
in step 201 , in which L*Pdiscrete channel estimations h 0 lo h LPm1 are developed. The channel estimate may be obtained 
using any conventional technique, e.g., by using correlators tuned to the pilot channel. Each of discrete channel esti- 
mations is taken with a time spacing of the chip duration divided by P, for CDMA, or symbol duration divided by P for 
TDMA. Also, in step 201 , the background noise plus interference power a 2 and the power of the downlink from the 
base station to the terminal a 2 are determined in the conventional manner. Fn step 203, the weights employed by joint 
equalizer 1 09 are determine^ as will be described fully further hereinbelow. 

[0022] Next, in step 205, a set of P samples from each antenna is obtained. Thereafter, in step 207, the determined 
weights of joint equalizer 1 09 are applied to the samples from each antenna by joint equalizer 109. The samples are 
then despread by despreader 111 , if CDMA was employed, in optional step 209 In optional step 211 , conventional soft 
mapping of the symbols to soft bits is performed, and the resulting soft bits are supplied as an output for use by a 
decoder. Thereafter, conditional branch point 213 tests to determine whether one coherence time has elapsed since 
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the previous execution of step 201. If the test result in step 213 is NO, indicating that the channel is believed to still 
remain substantially the same as when it was last estimated, control passes to step 205, and the process continues 
as described above. If the test result in step 213 is YES, indicating that sufficient time has passed such that the channel 
may have changed enough so as not to be considered substantially the same as when it was last estimated, control 
5 passes back to step 201 and the process continues as described above. 

[0023] FIG. 3 shows in more detail the process of step 203 by which the weights employed by joint equalizer 1 09. 
Note that the process requires the use of several matrices, and the dimensions of the various matrices are given in 
Table 2. 



10 



Table 2 



15 



20 



25 



Matrix 


Dimension 


r p 


EPXEP 


r pp 


NEP X NEP 




EP X (E+L-1) 


r(H) 


NEP X M(E+L-A) 


r(Hj 


NEPxM(E+L-1) 




(E+L-1 )X1 


A 


MX M(E+L-1) 




1xM(E+M) 


w 


MX NEP 




1 X NEP 


Q 


MX M 



[0024] In step 301 , the channel estimates h 0 to /7 LP . r for each transmit and receive pair obtained in step 201 are 
arranged in a respective matrix h n m as shown in equation 1 . In step 303, matrix convolution operator r(h n m ) is then 
formed for each respective matrix h nm , as shown in equation 2. Thereafter, in step 305, MIMO convolution operator 
r(H) is then formed from the various matrix convolution operators as shown in equation 3. 



30 



Km = 



35 



equation (1) 



\L-\)P 



40 



45 



r(h,„) = 



50 



0, h m 



equation (2) 
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w 



r(H) = 



r(h u ) r(h,,) - r(h, M ) 



r(h w ) 

equation (3) 



[0025] In step 307, delay vector e d is formed as shown In equation 4. Delay vector e d is a one dimensional vector 
with E+L-1 elements which are all zero except for the single value at the E+L-1 -d location, which has a value of 1 . A 
15 typical value for d, which is selectable by the implementor, is such that the location which has a value of 1 is at the 
center of the vector. The purpose of delay vector e d , is to impose the overall equalizer delay d onto the equalizer. , 
Thereafter, in step 309, delay matrix A is computed from equation 5, in which \ M is an identity matrix of size MxM. 



20 



e d = [0-10 0] 



equation (4) 



A=l w «e d 



equation (5) 



25 



[0026] Finally, in step 311, equalizer weight matrix, W, is computed in accordance with equation 6, in which X H means 
the Hermitian transpose of X, which is the complex conjugate transpose of the vector or matrix X. 



w = Ar(H) H (r(H) H r(H + r)R p / 



equation (6) 



30 



35 



[0027] In one embodiment of the invention, execution of step 205 is such that the P samples r of antenna n are 
initially arranged as a vector c shown in equation 7, where k\s the current received chip time index if CDMA is employed, 
or the symbol time index if CDMA is not employed. E consecutive in time vectors c for antenna n are then arranged 
as shown in equation 8, and the E consecutive in time vectors c for all of the N antennas are further arranged as shown 
in equation 9, forming a vector of received samples at time k. 



40 



45 



equation (7) 



50 



,(*) = 



55 



c„(*-£+l) 



equation (8) 
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10 



•(*)- 



'■(*) 



/» (*). 
equation (9) 



/5 



[0028] The application of the determined weights by joint equalizer 109 in step 207 may be performed as shown in 
equation 1 0, where y(/c) is a resulting vector of size M x 1 which contains the equalized chips if CDMA is employed, 
or symbols if CDMA is not employed. 



50 



y{k) = Wr(/c) 



equation (10) 



20 



25 



30 



35 



[0029] FIG. 4 shows another exemplary process, in flow chart form, for the overall operation of system of FIG. 1 . 
The version of the process shown in FIG. 4 is similar to that shown in FIG. 2, but it is for embodiments of the invention 
that include optional spatial whitening filter 113 (FIG. 1) and optional a posteriori probability (APP) metric processor 
115. Unless otherwise noted, all variables and parameters employed in the process of FIG. 4 are the same as described 
for FIG. 2. As with the process of FIG. 2, prior to performing the process of FIG. 4, initial values of parameters M, N, 
L, P, E, and c/must be determined. Additionally, prior to performing the process of FIG. 4 it is necessary to determine R pp . 
[0030] The process of FIG. 4 is executed periodically, with a periodicity that preferably does not exceed the coherence 
time, T co , which is the time duration for which the channel properties are substantially constant. The process is entered 
in step 401, in which L*P discrete channel estimations h Q \o /j lp _j are developed. Also in step 401, the background 
noise plus interference power o 2 and the power of the downlink from the base station to the terminal a 2 are determined 
in the conventional manner. In step 403, the weights W employed by joint equalizer 109 are determined, as described 
hereinabove. Additionally, in step 403, the effective channel matrix H eff and, optionally, spatial whitening filter Q are 
determined, in accordance with an aspect of the invention. More specifically, H eff is determined as shown in equation 
11, and Q is determined as shown in equation 12. 



H en (m,n) = e n T Wr(H)e m m = 1 :M, n = 1 :N 



equation (11) 



40 



45 



Q = (WW H )* 1/2 



equation (12) 



[0031] Next, in step 405 a set of P samples from each antenna is obtained. Thereafter, in step 407, the determined 
weights of joint equalizer 1 09 are applied to the samples from each antenna by joint equalizer 1 09. The samples are 
then despread by despreader 1 1 1 , if CDMA was employed, in optional step 409. In step 411, whitening filter Q is applied 
to the despread, equalizer outputs of step 409. In step 413, APP softbits are computed, in accordance with an aspect 
of the invention, using equations 1 3 and 1 4. The softbits are the output of the process of FIG. 4, and they may be made 
available to a decoder, e.g., a turbo decoder. 



55 
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5 




10 



equation 13 



15 




J 



equation 14 



20 



[0032] Thereafter, conditional branch point 415 tests to determine whether one coherence time has elapsed since 
the previous execution of step 401 . If the test result in step 415 is NO, indicating that the channel is believed to still 
remain substantially the same as when it was last estimated, control passes to step 405, and the process continues 
as described above. If the test result in step 41 5 is YES, indicating that sufficient time has passed such that the channel 
25 may have changed enough so as not to be considered substantially the same as when it was last estimated, control 
passes back to step 401 and the process continues as described above. 

[0033] FIG. 5 shows another exemplary embodiment of a multiple-input multiple-output (MIMO) wireless system in 
which the bit error rate floor caused by time dispersion is reduced by employing a joint minimum mean square error 
(MMSE) equalizer for all of the respective transmit antenna - receive antenna pairings that are possible in the MIMO 

30 system, in accordance with the principles of the invention. Shown in FIG. 5 are a) transmitter 501 ; b) transmit antennas 
503; including transmit antennas 503-1 through 503-M; c) receive antennas 505; including receive antennas 505-1 
through 505-N; d) receiver front-end processor 507; e) buffer-subtractor 521 ; f) joint equalizer 523; g) optional de- 
spreader 525; h) soft bit mapper 527; i) space-time regenerator 529; j) order controller 531 and k) switch 533. 
[0034] Transmitter 501 is a MIMO transmitter, e.g., one in which an original data stream is divided into substreams 

35 and each resulting substream is transmitted as a modulated radio signal via an individual one of transmit antennas 
503. The transmitted signals pass to the receiver over a time dispersive channel such that signals from each transmit 
antenna 503 reach each of receive antennas 505. 

[0035] Receive antennas 505 convert the radio signals impinging upon them into electrical signals, which are supplied 
to receiver front-end processor 507. Receiver front-end processor 507 operates conventionally to produce a stream 
40 of binary numbers representing samples of the radio signals received at antennas 505. Typically receiver front-end 
processor 507 performs radio frequency downconversion, filtering, sampling, and analog-to-digital conversion. The 
resulting samples are provided to buff er-substractor 521 . 

[0036] Buffer-subtractor 521 is shown in more detail in FIG. 6. Buffer subtractor 521 includes buffer 601 , memory 
603, and subtractor 605. Buffer 601 stores a time consecutive set of samples from each of antennas 505 as those 

^5 samples become available from front-end processor 507. Buffer 601 supplies the samples it receives from front-end 
processor 507 to memory 603 when an entire set is stored, i.e., the contents of buffer 601 are quickly dumped to 
memory 603. Once the set of samples are stored in memory 603 they may be independently accessed and sent to 
joint equalizer 523 (FIG. 5). Additionally, space time samples from a second input of buffer-subtractor 521 are supplied 
to subtractor 605. Subtractor 605 is cable of forming the difference between a specified location in memory 603 and 

so a second input to buffer-subtractor 521 . The resulting difference is stored in the specified location in memory 603. 
[0037] Joint equalizer 523 performs M passes through memory 603 with a different equalizer weight w m in each 
pass. Each weight is chosen to emphasize transmit antenna m and to suppress transmit antennas im- 1 through M. 
[0038] The output of joint equalizer 523 is supplied to conventional soft bit mapper 527, via optional conventional 
despreader 525 if CDMA is employed. The soft bits developed by soft bit mapper 527 are then supplied as an output, 

55 such as may be used by a decoder, e.g., the well known "Turbo decoder". 

[0039] The same output that is supplied to conventional soft bit mapper 527 may also be supplied to space-time 
regenerator 529 via switch 533. Space-time regenerator 529 forms a set of time consecutive samples for each receive 
antenna assuming the soft symbol is correct. In other words, assuming a particular soft symbol had been the actual 
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symbol transmitted by a particular transmit antenna, space-time regenerator 529 creates the corresponding effect that 
such a symbol would have caused on each of receive antennas 505 given the channel characteristics. Operation of 
space-time regenerator 529 will be explained more fully hereinbelow. 

[0040] Order controller 531 determines, based on channel estimates, the signal from which transmit antenna will be 

5 processed at any particular time, as will be explained more fully hereinbelow. 

[0041] FIG. 7 shows an exemplary process, in flow chart form, for the overall operation of system of FIG. 5 in which 
switch 533 is connected between despreader 525 and space-time regenerator 529. Prior to performing the process of 
FIG. 7, initial values of parameters M, N, L, P, E, and dmust be determined. What each of these parameters represents 
is listed in Table 1 . Additionally, prior to performing the process of FIG. 7 it is necessary to determine R pp , which is 

10 samples of the autocorrelation of the chip pulse shape r(t) t when CDMA is employed, or the symbol pulse shape 
autocorrelation, when CDMA is not employed. 

[0042] The process of FIG. 7 is executed periodically, with a periodicity that preferably does not exceed the coherence 
time, T co . The process is entered in step 701, in which L*P discrete channel estimations h 0 \o /? LP . f are developed. 
The channel estimate may be obtained using any conventional technique, e.g., by using correlators tuned to the pilot 

15 channel. Each discrete channel estimation is taken with a time spacing of the chip duration divided by P, for CDMA, 
or symbol duration divided by P for TDM A. Also in step 701 , the background noise plus interference power o 2 and the 
power of the downlink from the base station to the terminal o 2 are determined in the conventional manner. " 
[0043] In step 703, order controller 531 determines, according to equations 1 5 and 1 6, the order in which the signals 
from the various transmit antennas will be processed, with the signal from a respective transmit antenna being proc- 

20 essed for each execution of joint equalizer 523. Sort is a function that rearranges the elements of vector P so that they 
run from largest to smallest and orderls a list of all the antenna transmit antenna numbers as they should be processed 
by joint equalizer 523. It is preferable to process the so-called "strong" signals first. However, the particular character- 
istic, or set of characteristics, which are used to define the "strength" of a signal is at the discretion of the implementer. 
In the particular embodiment shown herein, estimated signal powers are employed as the strength. Further note that 

25 although herein the order for all the antennas is determined simultaneously, those of ordinary skill in the art will readily 
recognize that it is possible to successively determine which antenna to process. 



30 



P = \\diag(W r(H))f 
[P, order] = sort{P) 



equation (15) 
equation (16) 



35 



[0044] In step 705, the equalizer weights for the particular antenna currently being processed, m, as specified by 
the order, is determined according to equations 17 and 18, in which delay vector a m is the ni h row of delay matrix A 
which was described hereinabove. 



40 



w m = a m T(H m ) H |r(H m ) H r(H m ) + ^R 



2 PP 



equation (17) 



45 



50 



r(H.)« 



o r(h u ) o 
o r(h 2 ,) o 



.o r(iy 2 ) o 

equation (18) 



55 



[0045] Note that for each iteration m of equation 18, columns corresponding to order(1 ) through order(m-1 ) are set 
to the block zero matrix 0. Note that order is vector which contains a listing of the M antenna numbers in the order in 
which they will be processed. Order may be the result of the well known function sort of MatLab®. Doing so accounts 
for the fact that signals from transmit antennas 1 through m-1 have already been subtracted from the signal remaining 
to be processed for this set of samples. 
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[0046] Next, in step 707 a set of samples that span at least the duration of a data symbol is obtained from each 
receive antenna. Thereafter, in step 709, counter variable m is initialized to 1 . In step 71 1 , the weights determined for 
joint equalizer 523 in step 705 are applied to the samples from each receive antenna by joint equalizer 523, in accord- 
ance with equation 19. 

y(k) = w m T r(k) equation (1 9) 



[0047] The samples are then despread by despreader 525, if CDMA was employed, in optional step 71 3. In optional 
10 step 715, conventional soft mapping of the symbols to soft bits is performed, and the resulting soft bits are supplied 
as an output for use by a decoder. 

[0048] In step 71 7, samples are produced by space-time regenerator 529 according to equations 20 and 21 . 



15 



20 



s(kG) . ' 
s(kG+l) 

s(kG+G-l) 



equation (20) 



25 



30 



r(h^) 

equation (21) 



35 [0049] In step 719 the output of space-time regenerator 529 is subtracted from the contents of memory 603 (FIG. 
6), as shown by equation 22. 



rW = r(/c)-y m (*) equation (22) 

40 

[0050] Thereafter, conditional branch point 721 tests to determine if m is equal to M. If the test result in step 721 is 
NO, indicating that not all of the transmit antennas have yet had their signal contribution processed, m is incremented 
in step 723. Thereafter, control passes back to step 711 and the process continues as described above. If the test 
result in step 712 is YES, indicating that all of the transmit antennas have had their signal contribution processed, 

45 control passes to conditional branch point 725, which tests to determine whether one coherence time has elapsed 
since the previous execution of step 701 . If the test result in step 725 is NO, indicating that the channel is believed to 
still remain substantially the same as when it was last estimated, control passes to step 707, and the process continues 
as described above. If the test result in step 725 is YES, indicating that sufficient time has passed such that the channel 
may have changed enough so as not to be considered substantially the same as when it was last estimated, control 

50 passes back to step 701 and the process continues as described above. 

[0051] In an alternative configuration of FIG. 5, switch 533 is connected to the output of an error correction decoder 
which is either directly or indirectly connected to soft bit mapping 527. Doing so may improve performance. 
[0052] FIG. 8 shows a particular embodiment of joint equalizer 109, in which the equalization is performed in the 
discrete frequency domain, in accordance with an aspect of the invention. Shown in FIG. 8 making up joint equalizer 

55 109 are a) fast Fourier transform (FFT) processors 801 , including FFT processors 801 -1 through 801 -N, where N is 
the number of receive antennas; b) channel estimators 803, including channel estimators 803-1 through 803-N; c) fast 
Fourier transform (FFT) processors 805, including FFT processors 805-1 through 805-N*M, where M is the number of 
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transmit antennas; d) MMSE detection per frequency bin processor 807; and e) inverse fast Fourier transform (IFFT) 
processors 809, including IFFT processors 809-1 through 809-M. 

[0053] Each of FFT processors 801 receives from front-end 1 07 a signal of time domain digital samples correspond- 
ing to a respective one of receive antennas 105-N and performs the FFT algorithm on a set of consecutive samples 

5 to convert the time domain samples to samples in the discrete frequency domain, r n (<o), where a) is a particular discrete 
frequency. The number of samples F is at the discretion of the implementer based on a tradeoff between the perform- 
ance and complexity of FFT processors 801 . Typically the number of samples is a power of 2, e.g., 1 28, although the 
[jiiore samples the more accurate the equalization process will be. The possible values for co are determined as w = 
— , where n ranges from 0 to F-1 . 

10 ffie resulting discrete frequency samples for each receive antenna are supplied to MMSE detection per frequency bin 
processor 807. 

[0054] Each of channel estimators 803 also receives from front-end 107 a signal of time domain digital samples 
corresponding to a respective one of receive antennas 1 05-N and performs a channel estimate for the channel between 
its respective receive antenna and each of the M transmit antennas, thereby producing M channel estimates. Each 
15 channel estimate is a series of complex numbers that defines the impulse response of the channel. Fast Fourier trans- 
form (FFT) processors 805 each converts a respective channel estimate into the discrete frequency domain represen- 
tation thereof, and supplies the resulting discrete frequency domain representation of the channel estimates h n m (w) 
to MMSE detection per frequency bin processor 807. 

[0055] MMSE detection per frequency bin processor 807 performs the equalization in the frequency domain by com- 
20 puting 

z(o))=(H(co) H H(o>)+a 2 l)" 1 H(a>) H r(o>) equation (23) 

25 where 



Km W 



r x (a>) 




equation (25) 



30 



E(a>) = 



h u(*>) 



35 



equation (24) 



2 °n 

a equation (26) 

equation (26) 

55 

2 

a is the background noise plus interference power, 

c e 2 X is the power of the downlink signal from the base station to the terminal, and 
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I is the identity matrix. 

[0056] Each of the resulting M components of resulting vector z(a>) are then inverse frequency transformed from the 
discrete frequency domain into the time domain by inverse fast Fourier transform (IFFT) processors 809. The time 
domain equalized outputs are then supplied as the output of joint equalizer 1 09. 

[0057] FIG. 9 shows a particular embodiment of joint equalizer 109 in which the equalization is calculated in the 
discrete frequency domain and applied in the time domain, in accordance with an aspect of the invention. Shown in 
FIG. 9 making up joint equalizer 109 are a) matrix finite impulse response (FIR) filter 901 ; b) channel estimators 903, 
including channel estimators 903-1 through 903-N; c) fast Fourier transform (FFT) processors 905, including FFT 
processors 905-1 through 905-N*M, where M is the number of transmit antennas; d) MMSE tap weight calculator 907; 
and e) inverse fast Fourier transform (IFFT) processors 909, including IFFT processors 909-1 through 909- N*M 
[0058] Matrix finite impulse response (FIR) filter 901 continuously receives from front-end 1 07 a signal of time domain 
digital samples corresponding to a respective one of receive antennas 1 05-N. The number of taps of matrix FIR filter 
901 is at the discretion of the implementer based on a tradeoff between performance and complexity. Typically the 
number of samples is a power of 2, e.g., 128. 

[0059] Each of channel estimators 903 also receives from front-end 107 a signal of time domain digital samples 
corresponding to a respective one of receive antennas 1 05-N and performs a channel estimate for the channel between 
its respective receive antenna and each of the M transmit antennas, thereby producing M channel estimates. Each 
channel estimate is a series of complex numbers that defines the impulse response of the channel. Fast Fourier trans- 
form (FFT) processors 905 each converts a respective channel estimate into the discrete frequency domain, and sup- 
plies the resulting discrete frequency domain representation of the channel estimates h n m (a>) to MMSE tap weight 
calculator 907. The number of samples employed by FFT processors 905 for each conversion should be the same as 
the number of taps in matrix FIR filter 901 . 

[0060] MMSE tap weight calculator 907 develops frequency domain representations of the weights necessary to 
perform the equalization in the time domain by computing 



where H(o>), a 2 are defined as explained hereinabove in connection with FIG. 8. 
[0061] Each of the resulting M components of resulting vector S((o) are grouped by frequency into frequency vectors, 
and are then the frequency vectors are inverse frequency transformed from the discrete frequency domain to become 
filter weights in the time domain by inverse fast Fourier transform (IFFT) processors 909. The weights are then supplied 
to matrix FIR filter 901 which utilizes them to perform equalization in the time domain on the signals received from 
front-end 1 07 as shown by equation 28, 



where y(A) is the vector output at time k, y having M components — one for each transmit antenna — , Sy is the 
MxN filter matrix for delay j, which is the inverse Fourier transform of S(a)), r (/c) is the vector input signal which is 
received by matrix FIR filter 901 as defined in equation (9), and F is the number of samples taken for each FFT. 
[0062] Given the foregoing, those of ordinary skill in the art will readily recognize that other equalizer algorithms 
which approximate the operation and performance of MMSE, such as least mean square (LMS), recursive least squares 
(RLS), or minimum intersymbol interference (ISI) subject to an anchor condition, can be employed in a joint manner, 
e.g., in a space manner, in the implementation of joint equalizer 109. 

[0063] Those of ordinary skill in the art will readily recognize that the techniques of the instant invention may be 
employed in systems in which the various transmit antennas are transmitting at different data rates, e.g., using different 
encoding rates and/or transmit constellations, such as quaternary phase-shift keying (QPSK) or 16-ary quadrature 
amplitude modulation (16-QAM). In such a situation, if the embodiment of FIG. 5 is employed, order controller 531 is 
not required because the antennas must always be processed in an a priorily determined order. 



S(o))=(H(o)) H H(a))+a 2 l)' 1 H(o) 



H 



equation (27) 




equation (28) 
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Claims 

1 . A method for compensating for time dispersion in a receiver of a wireless system that has a plurality of transmit 
antennas and a plurality of receive antennas, the method comprising the steps of: 

5 

receiving samples for each receive antenna; 

determining a joint equalizer solution using channel information for at least one pairing of at least one of said 
transmit antennas and said receive antennas; and 

applying said determined joint equalizer solution to said received samples from at least one of said receive 
10 antennas to develop equalized samples. 

2. The invention as defined in claim 1 wherein said joint equalizer solution is a joint minimum mean square error 
(MMSE) solution. 

15 3. The invention as defined in claim 1 further comprising the step of estimating a channel for at least one pairing of 
at least one of said transmit antennas and said receive antennas. 

4. The invention as defined in claim 1 wherein at least two of said transmit antennas transmit antennas transmit at 
least one of different rates; and 

20 different transmit constellations. 

5. The invention as defined in claim 1 further comprising the step of performing soft bit mapping using a version of 
said equalized samples. 

25 6. The invention as defined in claim 5 wherein said version of said equalized samples are despread samples. 

7. The invention as defined in claim 5 wherein said step of performing soft mapping further comprises the step of 
spatial whitening said version of said equalized samples. 

30 8. The invention as defined in claim 5 wherein said step of performing soft mapping further comprises the step of 
performing a posteriori probability (APP) metric processing on said version of said equalized samples. 

9. The invention as defined in claim 1 wherein said determining and applying steps are iterated multiple times over 
a symbol period, one iteration for each one of said transmit antennas, and said method further comprises, for each 

35 iteration, the steps of: 

generating a representation of signals that would have arrived had a particular symbol for a currently being 
processed transmit antenna had been transmitted; and 

subtracting said representation from said samples received for each receive antenna. 

40 

10. A receiver for use in a multiple-input multiple-output (MEMO) system in which a plurality of signal detectors receive 
signals transmitted by a plurality of signal sources, said receiver comprising: 

a joint equalizer that develops a joint equalizer solution using channel information for at least one pairing of 
45 said at least one of said signal sources and said signal detectors and supplies as an output a signal that 

includes at least said equalizer solution applied to a signal received by at least one of said signal detectors; and 
a soft bit mapper for developing soft bits from said joint equalizer output. 

11. The invention as defined in claim 10 wherein said joint equalizer solution is a joint minimum mean square error 
50 (MMSE) equalizer solution. 

12. The invention as defined in claim 11 wherein said joint equalizer further comprises: 

a first plurality of fast Fourier transform processors, each of said fast Fourier transform processors being 
55 coupled to receive a respective input corresponding to a signal received by one of said signal detectors and 

supplying as an output a discrete frequency domain representation thereof; 

a plurality of channel estimation units each of which is coupled to receive a respective input corresponding to 
a signal received by one of said signal detectors which develops a channel estimate for each channel between 
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10 



each respective signa! source and each respective signal detector; 

a second plurality of fast Fourier transform processors, each of said second plurality of fast Fourier transform 
processors coupled to receive a respective input corresponding to a channel estimate for a respective one of 
said channels between said signal sources and said signal detectors and supplying as an output a discrete 
frequency domain representation thereof; 

an MMSE detection per frequency bin processor coupled to receive as inputs said outputs from said first 
plurality of fast Fouriertransform processors and from said second plurality of fast Fouriertransform processors 
to produce a discrete frequency domain representation of an application of a joint minimum mean square error 
(MMSE) equalizer solution to said signals received by each of said signal detectors; and 
a plurality of inverse fast Fourier transform processors which convert said representation of an application of 
a joint minimum mean square error (MMSE) equalizer solution to the time domain. 



13. The invention as defined in claim 12 wherein M is the number of signal sources and N is the number of signal 
detectors, and 

15 wherein said MMSE detection per frequency processor performs the equalization in the frequency domain 

by computing 



z(a>) = (H(o)) H H(a>)+a 2 l)* 1 H(G>) H r(o) 



20 



where 



25 



30 



KM 



KM 



>>w(®) 



Km M. 



35 



40 



■(•)- 



r x {a>) 
r 2 (a>) 



45 



50 



2 

2 °n 
° ="2 
*x 

a 2 is the background noise plus interference power, 

0 0 is the sum of the power received by all said signal detectors from all of said signal sources, 
each r(w) is said output of a one of said first plurality of fast Fourier transform processors, 
each h((o) is said output of a one of said second plurality of fast Fouriertransform processors, 

1 is the identity matrix, and 

X H means the Hermitian transpose of X, which is the complex conjugate transpose of the vector or matrix X. 



14. The invention as defined in claim 11 wherein said joint equalizer further comprises: 



55 



a plurality of channel estimation units each of which coupled to receive a respective input corresponding to a 
signal received by one of said signal detectors which develops a channel estimate for each channel between 
each respective signal source and each respective signal detector; 

a plurality of fast Fouriertransform processors, each of said plurality of fast Fouriertransform processors being 
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coupled to receive a respective input corresponding to a channel estimate for a respective one of said channels 
between said signal sources and said signal detectors and supplying as an output a discrete frequency domain 
representation thereof; 

an MMSE tap weight calculator coupled to receive as inputs said outputs from said plurality of fast Fourier 

transform processors to produce a discrete frequency domain representation of a joint minimum mean square 

error (MMSE) equalizer solution to said signals received by each of said signal detectors; 

a plurality of inverse fast Fourier transform processors which convert said representation of a joint minimum 

mean square error (MMSE) equalizer solution to matrices of filter coefficients in the time domain; and 

a matrix finite impulse response (FIR) filter coupled to apply said matrices of filter coefficients in the time 

domain to said signals received by said signal detectors. 



15. The invention as defined in claim 14 wherein M is the number of signal sources and N is the number of signal 
detectors, and wherein said MMSE equalizer solution is developed by computing: 



15 



S(o))=(H((i)) H H(o))+a 2 l)" 1 H(a)) H 



where 



20 



25 



h u(*>) 



30 




a is the background noise plus interference power, 
35 a 0 is the sum of the power received by all said signal detectors from all of said signal sources, 

each h(o>) is said output of a one of said plurality of fast Fourier transform processors, 
I is the identity matrix, and 

X H means the Hermitian transpose of X, which is the complex conjugate transpose of the vector or matrix X. 

40 16. The invention as defined in claim 15 wherein said matrix FIR filter applies said matrices of filter coefficients in the 
time domain to said signals received by said signal detectors by computing 



where y(A) is the vector output at time k } y having one component for each of said signal sources, 
Sy is a MxN filter matrix for delay j, which is the inverse Fourier transform of S(o>), 
so r(A) is a vector of received samples at time k, and 

F is the number of samples taken for each fast Fourier transform by each of said plurality of fast Fourier 
transform processors. 

17. The invention as defined in claim 10 wherein said joint equalizer develops said joint equalizer solution as a function 
55 of estimates of the channels between each of said signal sources and said signal detectors. 

18. The invention as defined in claim 11 wherein said joint equalizer develops said joint minimum mean square error 
(MMSE) equalizer solution by computing 
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W = A 



J o 2 V 1 

r(H) w (r(H)1r(H) +1 n R pp 



5 where 

r(H) is a MIMO convolution operator, 

X H means the Hermitian transpose of X, which is the complex conjugate transpose of the vector or matrix X, 
A is a delay matrix, 

a 2 is the background noise plus interference power, 
10 a 0 is the sum of the power received by all said signal detectors from all of said signal sources, 

noise covariance R pp , 

and said joint equalizer applies said joint minimum mean square error (MMSE) equalizer solution by com- 
puting 

y(A) = Wr(/r) 

15 where r(kj is a vector of received samples at time fr. 

19. The invention as defined in claim 10 wherein said joint equalizer solution is one from the group consisting of: a 
joint least mean square (LMS) solution, a joint recursive least squares (RLS) solution, or a joint minimum inter- 
symbol interference (ISI) subject to an anchor condition solution. 

20 
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FIG. 3 
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FIG. 4 
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FIG. 6 
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